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Abstract. 

Assuming that general relativity is the correct theory of gravity in the strong field 
limit, can gravitational wave observations distinguish between black hole and other 
compact object sources? Alternatively, can gravitational wave observations provide a 
test of one of the fundamental predictions of general relativity? Here we describe a 
definitive test of the hypothesis that observations of damped, sinusoidal gravitational 
waves originated from a black hole or, alternatively, that nature respects the general 
relativistic no-hair theorem. For astrophysical black holes, which have a negligible 
charge-to-mass ratio, the black hole quasi-normal mode spectrum is characterized 
entirely by the black hole mass and angular momentum and is unique to black holes. 
In a different theory of gravity, or if the observed radiation arises from a different 
source (e.g., a neutron star, strange matter or boson star), the spectrum will be 
inconsistent with that predicted for general relativistic black holes. We give a statistical 
characterization of the consistency between the noisy observation and the theoretical 
predictions of general relativity, together with a numerical example. 
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1. Introduction 

The formation of a black hole is the ultimate manifestation of strong field gravity. 
During the late stages in the aspherical formation of an astrophysical black hole the 
gravitational waves emitted are dominated by a set of quasi-normal modes (QNMs) 
[H El 0]: exponentially damped sinusoids, whose frequency and damping times are 
characteristic of the black hole's mass and angular momentum (And electric charge, 
as well; however, astrophysical black holes, which are our interest here, have negligible 
charge-to- mass ratio). Should we observe a QNM from a black hole and know also which 
particular normal mode we are observing we can determine, from the mode's frequency 
and damping time, the black hole mass and angular momentum [HE]. 

If, on the other hand, we observe several different QNMs from the same source 
and find that they are inconsistent with the spectrum predicted by general relativity 
in the sense that they cannot be explained by a single value of the mass and angular 
momentum we may infer that we are not observing a black hole. Here we develop this 
observation into a experimental test of the existence of black holes or, alternatively, 
general relativity. 

Future observations by the Laser Interferometer Space Antenna (LISA) [I] offer us 
a different perspective on this question. LISA is expected to observe mergers of compact 
objects with masses in the range 10 6 — 1O 8 M 0. In our present understanding, these 
compact objects can only be black holes. Observations by LISA of QNMs inconsistent 
with black holes would also be a test of the general relativistic no-hair theorem since an 
inconsistency in this mass range with black hole sources would indicate that physical 
scales other than mass and angular momentum were involved in the generation of the 
radiation. 

In either sense the test described here is of general relativity based on gravitational 
wave observations. Eardley et al. (9^ proposed the first test of general relativity 
using gravitational wave observations. They investigated the polarization modes of 
gravitational waves in various metric theories of gravity and described how to identify 
the modes experimentally and use those observations to identify the spin content of 
dynamical gravity. The first actual test of general relativity relying on its prediction 
of the existence of gravitational waves was made by Taylor and Weisberg [HI]. They 
described how the observed orbit and orbit decay of the Hulse- Taylor binary pulsar 
PSR B1913+16 led to a strong consistency check on the predictions of general relativity. 
Finn (TTj proposed a different test of the spin content of dynamical gravity, based on 
the possibility of a space-based detector in circumsolar orbit observing the induction- 
zone field associated with solar oscillations. Ryan ^2] has outlined how observations of 
the gravitational radiation from capture orbits of solar mass compact bodies about a 
supermassive black hole may allow the determination of certain multipole moments of 
the central hole, thereby testing the prediction of general relativity. More recently, Will 
[T3] and Finn and Sutton jUj IT5] have described tests of general relativity that bound 
the mass of the graviton, and Scharre and Will JH] and Fairhurst et al. 17 j have shown 
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how gravitational wave observations of pulsars may be used to bound the value of the 
Brans-Dicke coupling constant. 

The preceding tests can be grouped into three different classes. One set of tests, 
which includes |lU|IT^llTH ll6j. is based on energy conservation arguments: the observed 
evolution of a system or of the radiation from a system, is related to the energy loss 
expected owing to the radiation. A second class of tests, which includes jHJ HU ITTj . 
focuses on the observed polarization modes of the field. A third class, which includes 
[T3] , involves the frequency- dependent dispersion relationship associated with a massive 
graviton. The test described here is of a new class, based on the unique character of the 
radiation spectrum associated with a disturbed black hole. 

QNMs appear as solutions to the equations describing perturbations of a stationary 
black hole spacetime, subject to the boundary conditions of no in-going radiation from 
infinity and no up-coming radiation from the horizon. The perturbation equations 
describing Schwarzschild black holes were first described by Regge and Wheeler [TH] 
and Zerilli jTHlEDl- The first QNM solutions to these equations were found by Press 
[21]. Teukolsky found the corresponding perturbation equations for Kerr black holes 
[221 EH], and, with Press, first investigated their QNM solutions [21] . 

Damped sinusoidal motion is ubiquitous for systems approaching equilibrium and 
one expects that collapse or coalescence will lead, in any theory of gravity, to some form 
of QNM ringing. If we observe a QNM spectrum that is inconsistent with an isolated 
black hole, then there are two possibilities. On the one hand, general relativity may yet 
be correct, but we are not observing an isolated black hole approaching equilibrium" 1 ". 
Similarly, we may be observing the radiation arising from a compact body that is not 
a black hole — e.g., a neutron star [2E], a boson star j2£] or strange matter star |27j . 
whose QNM spectrum will be determined by the properties and configuration of the 
appropriate matter fields — or a black hole carrying a previously unknown macroscopic 
charge e.g., a dilaton field [2H]- On the other hand, general relativity may not be the 
correct theory of gravity in the strong field limit. Thus, while no single observation may 
rule out general relativity, a set of observations, each of a different source, none of which 
is consistent with an isolated black hole, could suggest the need to consider alternative 
theories of gravity in the strong- field limit. 

This paper is organized as follows. In section [2] we briefly describe the QNMs 
of a Kerr black hole and explain how the idealized observation of two or more modes 
in the absence of noise enables us to extract the mass and angular momentum of the 
black hole. Real boats, of course, rock, and section [3] generalizes the discussion to 
include experimental errors and describes how one can use noisy gravitational wave 
observations of QNMs to test general relativity. In section 0] we demonstrate, in a 
model numerical simulation, the use of this method as applied to LISA observations. 
Section E3 investigates the range to which we can expect LISA to observe sources strong 

+ In fact, the uniqueness theorems have only been proved for vacuum spacetimes and they are not true 
in the presence of arbitrary matter fields or radiation. Nevertheless, it would be a great surprise if the 
spacetime in the vicinity of a black hole is not close to Kerr in some approximate sense. 
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enough for this test to be applied. We conclude in section |U] with a summary of our 
main results. 

2. Ideal observations 

2.1. quasi-normal modes of Kerr black holes 

Following the aspherical collapse to a black hole, one expects that the final black hole 
can be described as the perturbation of a stationary Kerr hole. The dominant part of the 
gravitational waves emitted as the black hole settles down can be described as a sum over 
a countably infinite set of damped sinusoids, each characterized by an amplitude, phase, 
frequency and damping time. (At still later times, the radiation will be dominated by 
power-law tails arising from the backscatter of radiation off the spacetime curvature in 
the neighborhood of the black hole PI El El; however, here we are interested in the earlier, 
and higher amplitude, QNM ringing.) In this sub-section we review those properties of 
the black hole QNM spectrum that are important for our investigation; more detailed 
examinations of the spectrum itself can be found in j2Hl EDI Ell E2] ■ 

Gravitational wave detectors respond to a linear combination of the radiation in the 
two polarization modes of the incident gravitational waves. In terms of the transvere 
traceless (TT) [33] gauge metric perturbation the observable quantity h(t) may be 
written, for QNMs, in the form 



where the summation indices characterize the particular mode, which is related to the 
angular dependence of the mode amplitude and phase on a sphere of constant (Boyer- 
Lindquist) radius about the black hole through I and m, and the "harmonic" through 
the index n: i = 2, 3 . . ., \m\ < £ and n = 1, 2 . . .. For the Schwarzschild geometry the 
symmetry is spherical, the appropriate decomposition of the metric perturbation is given 
by the usual spherical harmonics, and modes differing only in m are degenerate. For Kerr 
the symmetry is axisymmetric and the orthonormal decomposition of the perturbation 
is by spheroidal harmonics [23]. The amplitudes A nim and phases <t> n tm depend on the 
initial conditions and the relative orientation of the detector and the source; however, 
the complex frequency uj n g m depends only on the intrinsic parameters of the underlying 
black hole: i.e., its mass M and angular momentum aM 2 . (We assume that the black 
hole carries no significant electric charge.) 

For fixed a the complex frequency u n £ m scales as M _1 ; thus, we define the 
dimensionless frequency fl n £m, 




(1) 



? 




■m 



(2) 
(3) 
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where F nem and T n £ m are the real dimensionless frequency and damping time of the 
modes. The corresponding physical frequency f n e m and damping time r n g m are given by 

LOnim = 2nf n £ m + i/T n £ m = 2lT F nirn /M + ij (MT n i m ). (4) 

(We use geometrical units with G = 1 and c = 1.) The dimensionless VL n £ m (or F n £ m and 
T n £ m ) depend only on the also dimensionless black hole angular momentum parameter 
a. Figure H shows fl n £ m as a function of a for n — 1, 2 and £ = 2, 3, and |m| < £. 

From quasi-normal modes to testing relativity 

If we observe only one mode, characterized by its complex frequency u (cf. eq. HJ), what 
can we say about the underlying black hole? 

Corresponding to the observed u is the line Q = Mu, M e M>o, in the dimensionless 
Q plane (cf. eq.EJ). Such a line is shown in figure ^ This line will intersect some subset of 
the family of fl n e m curves, characteristic of black hole normal modes. Each intersection 
corresponds to a black hole mass M, angular momentum parameter a, and mode nlm 
consistent with the observed u. Knowing only / and r, then, we cannot uniquely identify 
the black hole mass and angular momentum, but we can reduce the possibilities to a 
(possibly countably infinite) set of (a, M) pairs. If we knew n£m as well, we would know 
a and M exactly. 

Now suppose that we observe two modes from the same black hole, each 
characterized by its own frequency and damping time. Figure shows, in schematic 
form, the line Mu for each of the two modes (denoted + and x ) and their intersection 
with several different Q n e m curves in the complex Q plane. Corresponding to each mode 
is a set of candidate (a, M) pairs that may describe the underlying black hole. Each 
candidate mass and angular momentum parameter is a point in the (a, M) plane, as 
shown in figure Eb- With two or more modes, there must be at least one common 
candidate mass and angular momentum. 

This is, in essence, our proposal for a test: interpreting the observation of several 
normal modes Uk, k > 2, as arriving from a single, general relativistic black hole, we 
conclude that general relativity is self-consistent if the observed Uk are consistent with 
at least one black hole (a,M). If no such (a, M) exists for the observed either we 
have observed something other than an isolated black hole or we have a contradiction 
with the predictions of the theory. 

(As an aside, it is possible (though unlikely) that we get more than one value 
of (a, M) consistent with the observed frequencies. This can happen if we have two 
mode pairs (ni,^,m 1 ;n 2 ,4) m 2) and (fix, £x, rhx, n 2 , £2, ^2) which give rise to the same 
frequency lu. In this case the observations would still be consistent with general relativity 
though we could not use that observation to measure M and a. The important point of 
our test is the existence of at least one (a, M) pair consistent with the observations.) 

Noise and other experimental realities ensure that there will be no exact agreement 
between the observed ujj, and a general relativistic black hole even if general relativity 
is correct. The challenge, then, in developing a practical test is to determine when the 



Black Hole Spectroscopy 



6 



Re(Q) 



Im(O) -0.3 - 



-0.4 - 



-G.6 



-0.2 - 



-0.5 - 



-0.1 - 








Figure 1. The dimensionless, complex QNM frequencies flnim for rotating, uncharged 
black holes. Each family of curves corresponds to one n£ pair, and each branch to a 
possible value of m. The large black dot at the base of each family is the Schwarzschild 
(a = 0) limit, where the frequencies are degenerate in m. This degeneracy is broken 
for a ^ 0, and the curves emanating from the dots give the QNM frequencies for Kerr 
black holes as a function of positive a for different m. In this figure a ranges from to 
0.9958, with the small diamonds marking the QNM frequencies for a = 0.4, 0.6, 0.8, 0.9, 
and 0.98. On this figure, an observation, corresponding to a (complex) frequency u>, 
is represented by the line = Mlu, parameterized by the (unknown) black hole mass 
M. Each intersection of this line with a QNM curve in dimensionless f2 represents a 
candidate n£m, M and a for the mode. 



differences between the candidate (a, M) pairs associated with the different observed 
modes are so great as to be statistically inconsistent with general relativity. In the next 
section we face this challenge. 

3. A test of relativity 

3.1. A reformulation of the test 

Before we discuss the role that noise plays in our analysis it is helpful to reformulate 
the test described in section EIH and figure 121 Consider an ordered iV-tuple of QNMs, 



Each Q may be regarded as a function that maps M and a to a set of observable 
frequencies 




(5) 




(6) 
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Figure 2. Here we show, in schematic form, several Vt n £ m {a) curves and their 
intersection with the lines Muji, M > 0, i = 1, 2, corresponding to two observed modes. 
These two lines we denote by + and x respectively, (b) The candidate (a,M) pairs 
determined in figure (a) are plotted here in the (a, M)-plane. The pairs belonging to 
u>i are denoted by +, the ones belonging to 0J2 by x. There is only one candidate 
(a, M) consistent with both observations — indicated by the overlapping + and x - 
and this is the actual mass and angular momentum of the underlying black hole. 



Each iV-tuple Q thus describes a two dimensional surface in the (2N + 2) -dimensional 
space S, 

S := (a,M,a>i, . . .,u N ), (7) 

with different A-tuples corresponding to different sets of N modes. (In section EHU we 
will understand the uj^ to represent observed QNM frequencies and damping times.) 
An observation oj consists of an A-tuple 

uj := (u 1 ,...,u N ). (8) 

The observation u) also corresponds to a surface in S. The observation is consistent with 
a black hole if the surface of constant uj intersects one of the surfaces Q. Figure El shows 
a low-dimensional projection of such an observation uj together with several surfaces 
(which appear as curves) for different iV-tuples Q. A moment's consideration should 
convince one that this new criterion is equivalent to the criterion formulated above in 
section 12.21 

In practice the situation is less ideal: noise distorts our observation, so that — even 
if we are observing black hole QNMs — the measured uj will not intersect a curve Q. 
In the remainder of this section we describe how our ideal test is made practical and 
meaningful for real observations. 

3.2. Confidence intervals and testing general relativity 

In a frequentist analysis, the observation, the sampling distribution, an ordering 
principle and a probability combine to determine a confidence interval. In this section 
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(a,Af)' 




Figure 3. A reformulation of the consistency criterion. A set of quasi-normal 
modes Q = {(nfe,ife,mfe) : k = 1, ... ,JV} corresponds to a surface in the (2iV + 2)- 
dimensional space depicted in this figure. A measurement w = (wi, . . . , ujn) is 
consistent with general relativity if the constant surface that is obtained by ranging 
over all (a, M) while keeping the frequencies uj fixed intersects at least one of the 
surfaces corresponding to one of the sets Q. This intersection is indicated in this figure 
by a dot. 



we use this construction to form a confidence region in the (a, M)-plane, given a noisy 
observation u>. 

We begin by reviewing the construction of a classical confidence interval for the 
one-dimensional case following jMj (alternatively, see e.g. We suppose that we 

make measurements of a random variable x from which a quantity fi is determined. 
The sampling distribution P(x\n) is the probability of making the observation x given 
a particular fi. Formally, an ordering principle is a function R(x\fi), which we use to 
identify a sub-interval J of a: according to 



The parameter r is chosen such that the region J(/x|r) encloses a fixed probability p: 



Given an observation xq, the probability-p confidence interval is the range of \i for which 
J{[i\r{p)) includes x as shown in figure EJ In an actual experiment, the choice of the 
value of the parameter p is made by the experimentalist. Typical choices are 90%, 95% 
or 99%. 

The choice of the ordering principle R(x\n) is a key ingredient in the construction 
of confidence intervals. Different choices will lead to different confidence intervals for the 
same observation: for example, one choice of ordering principle will always determine 
intervals of the form (— oo,x), while another choice will always determine intervals 
of the form (x, oo). Neither choice is a priori right or wrong. Here we will choose 
= P(x\fi) so that the intervals are given by level surfaces of the distribution 




(9) 




(10) 
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Figure 4. The construction of classical confidence intervals. A sampling distribution 
P(x\n), an ordering principle, and a probability p are needed to construct a confidence 
interval. The ordering principle is used to find the intervals J(/i) such that 
Ij(p) ^ x ^( x \^) = P- The classical confidence interval TZ is then given by the set 
of for which J(p) contains the measured value xq. 



P(x\fj). The main advantage of this ordering principle is that it is simple and it works 
in any dimension. Consider, for example, a two-component observation depending on 
one parameter a. There is, as before, a sampling distribution P(x, y\a) and an ordering 
principle R(x, y\a) = P(x,y\a). Confidence intervals can be defined in the same way as 
in the case of a one-dimensional observation; the interval J is now a two-dimensional 
region. Since this system is over-determined — we are now trying to determine one 
parameter a by measuring two quantities x and y — the measured x and y will have to 
satisfy additional constraints in order to give a non-vanishing confidence region. This 
is in fact precisely what happens in the black hole quasi-normal mode problem: any 
single measurement of uj can be explained by some (a, M), but a measurement of two 
or more uj can be simultaneously consistent with at least one (a, M) pair only if the 
no-hair theorem is true and the modes arise from a single black hole. 

We can now describe our test of relativity. Note that not all observations u> will lead 
to a non-empty confidence interval: i.e., for some u> there will be no (a, M) consistent 
with the observation. We frame our test in terms of the confidence interval we construct: 
if we make an observation uj for which there is no probability p confidence interval, then 
we say that the observed normal modes are inconsistent with an isolated black hole with 
confidence p. Conversely, if there does exist a probability p confidence interval, then we 
have verified that general relativity is self-consistent at this confidence level. 

Finally, we should point out an aesthetic flaw of our choice of ordering principle. 
The function P{x\p) is a density and, therefore, not invariant under a reparameterization 
of x. If we were to use a new parameter x' = f(x) for some smooth monotonic function 
/, the confidence region obtained for \i using a measurement of x may not coincide with 
the region obtained using a measurement of x' . In the one-dimensional case, there exists 
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another ordering principle based on the likelihood ratio (35] which is reparameterization 
invariant; however, we have not been able to generalize this to higher dimensions. While 
aesthetically displeasing, there is nothing wrong with the choice we have made, which 
is natural given the physical association of the parameters M and a with the black hole 
mass and angular momentum. 



3.3. Generalization to Quasi-normal Modes 

The generalization to QNM observations is straightforward. Each observation consists 
of N complex QNM frequencies u>k and associated amplitude signal-to-noise ratios 
Pk, which characterize both the amplitude of the signal at that frequency and the 
uncertainty in the determination of u>k (cf. [8J): 

uj := (ui, . . .,u N ) (11) 

P ■= (pi,---,Pn)- (12) 

For defmiteness suppose that u)k and pk are identified via maximum likelihood techniques 
jH]. There is a minimum signal-to-noise associated with each mode, which is set by the 
requirement that the observation must identify N modes. 

Observations u> corresponding to a black hole characterized by (a, M) and signal- 
to-noise pk are distributed according to the sampling distribution 

^ Probability of making observation ^ 
P(uj\a,M, Q, p) := u> given the actual iV-tuple Q . (13) 

^ and signals-to-noise p. j 

In general the sampling distribution depends upon the nature of the detector noise and 
the analysis procedure that identifies the modes Uk- For large signal-to-noise ratios it 
will generally reduce to a multivariate Gaussian in Ke(uk) and lm.(uk) and for smaller 
signal-to-noise ratios it can be determined via simulation. 

Now consider the region of the space S (cf. section 13. 1ft defined by 

P{u;\a,M,Q,p) > Po (14) 

with po such that 

P= [ P(u\a,M,Q,p)d 2N LU (15) 

JP(u>\a,M,Q,P)>PO 

for a fixed p. We say that the observation uj is consistent with a black hole if the 
actual observation uj is included in this region for some (a, M). Figure El illustrates the 
comparison of an observation with the region defined by equations ()14H15j) . It remains 
only to specify p. 

To help in specifying p it is useful to examine more closely its meaning. Suppose we 
have chosen a value of p. That value of p determines a confidence region. Now consider 
an ensemble of identical detectors, each observing simultaneously the same black hole 
event and its corresponding QNMs. The fraction of these observations that does not 
intersect the confidence region is the false alarm probability a(p), so-called because it 
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Figure 5. The construction of classical confidence intervals generalized to higher 
dimensions. Given a sampling distribution P, an ordering principle, and a probability 
p one can construct classical confidence regions R just as in the one-dimensional 
case. The difference here is that we are now trying to determine a small number of 
parameters (a, M) from a larger number of observations u> — (wi, . . . , wjv)- There are 
thus additional consistency conditions that need to be satisfied to obtain a non-empty 
confidence region R. This is the basis of our proposed test. 



is the probability that an observation will be falsely deemed to be inconsistent with a 
black hole. The probability a is a monotonic function of p; therefore, we can specify 
a in lieu of p. For observations whose characteristic frequency corresponds to masses 
greater than neutron star masses, which we are confident originate with black holes, we 
propose setting p so that a(p) — the probability of falsely rejecting the hypothesis that 
we have in fact observed a black hole — is small (e.g., less than 1%). In other words, 
the standard of evidence for declaring that we have discovered "new physics" should be 
high. 

The false alarm probability function a{p) will depend on the signal strength, as 
characterized by the signal-to-noise ratios; consequently, it will need to be determined 
on an observation- by-observation basis. Thus the calculation of a(p) by a Monte Carlo 
simulation is the final ingredient we need. In the next section we demonstrate the test 
through a numerical example where we calculate a(p). 
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4. A numerical example 

In the previous section we described a general procedure for testing general relativity by 
observing QNMs. In this section we explore its effectiveness numerically through a set 
of simulated observations drawn from a hypothetical black hole population population 
inspired by potential LISA observations, and a hypothetical population of non-black hole 
compact object sources, or NBHs. (We say "inspired" because, in fact, for the purpose 
of this analysis the observations are characterized entirely by the dimensionless signal- 
to-noise ratio and mode quality factor, with the dimensioned mode frequency simply 
setting a scale. Thus, the conclusions we reach are as valid for LISA observations as 
they are for observations at the same signal-to- noise with ground-based detectors.) 

For the BH observations we find the relationship between the false alarm probability 
a and the probability p that appears in equation (|15jl. For the NBH observations there 
are no "false alarms": every observation is of something not a black hole. Instead, there 
are false dismissals: observations that we mistakenly classify as consistent with a black 
hole. The probability of a false dismissal, denoted /3, depends on the choice of p or, 
alternatively, the choice of false alarm probability a(p) that we make for the purpose of 
defining the test. (The false dismissal probability depends also on how the spectrum of 
BHs and NBHs differ.) The smaller the false dismissal probability the more sensitive 
the test is to discovering "new physics" or identifying non-black hole sources. For the 
NBH observations we evaluate the false dismissal probability as a function of the false 
alarm probability. 



4-1. Simulating black hole QNM observations 

For defmiteness we focus on observations of two QNMs. For the purpose of 
illustration we consider black hole masses and angular momenta consistent with 
potential observations by the LISA detector jHJ. We first draw an (a, M) pair from 
the distribution 

P(a,M) = P(a)P{M) (16) 

1 for a e [0,0.986) 
' otherwise v ' 



P{M) oc 



M" 1 for M E (2.5 x 1O 5 M ,4.5 x 1O 8 M ) 
otherwise 



The range of M is determined by the frequency band where LISA is expected to be most 
sensitive; the range of a is determined by the maximum angular momentum expected 
of a black hole spun- up by thin-disk accretion [Sj. 

Corresponding to each (a, M) pair we choose the QNMs corresponding to (n = 
1,£ = 2,m = 2) and (n = 1,£ = 4, m — 4). We assign each mode the same signal- 
to-noise ratio, which we treat here as sufficiently large that the errors associated with 
the measurements are normally distributed with covariance matrix Cy equal to the 
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inverse of the Fisher information matrix 1^ (see e.g. [36 J as given in (HJ equation 
4.14]. This is in fact a mathematical lower bound — the Cramer-Rao bound — on the 
covariance matrix. We draw from this error distribution errors in the frequencies and 
damping times that we add to the "real" frequencies and damping times to determine 
the simulated observations: noisy QNM frequencies and damping times. 

Given this pair of QNM frequencies and damping times with errors we ask whether 
the two modes are in fact observationally distinguishable: if the frequencies and damping 
times are not sufficiently different then no real observation would ever result in the given 
pair. For instance, the five (n = 1, 1 = 2) modes are degenerate at a = 0; consequently, 
no matter how large the signal-to-noise ratio, if a is sufficiently small it is impossible to 
resolve these five modes observationally. 

To decide whether the two modes we are investigating are observationally 
distinguishable we invoke a "resolvability criterion" : denoting the frequencies (damping 
times) of the two modes as fx, f 2 (tx, r 2 ) we say that the two modes are distinguishable 



We discard any mode pair that does not satisfy this criterion. 

The result is an observation, which consists of a pair of signal-to-noise ratios and 
associated frequencies and damping times. (The observation does not include knowledge 
of black hole mass or angular momentum, or the n£m associated with the frequencies 
or damping times.) 

4-2. False alarm probability a 

For each simulated observation uj, constructed as described in section 14.11 we evaluate 
the smallest probability p = p min such that equations (fT3)l and (JT5Jl describe a region S 
that covers u) for some (a, M). The false dismissal fraction a(p) is the fraction of p m j n 
determinations that are greater than p: i.e., the fraction of BH observations that we 
would reject as originating from a black hole for threshold p. 

Ideally, in evaluating p we would consider every possible nim for each ujk- In 
practice, we consider only a finite subset of low-order (in both n and £) modes, 
corresponding to our expectation that these are the modes most likely to be excited 
to large amplitude. In our simulations we considered only modes corresponding to 
( n = 1J = 2, to = 0), (n = 1,£ = 2, to = 2), (n = 1,£ = 3, to = 3) and 
(n = 1,£ — 4, to = 4). Since for these simulations we observed two distinguishable 
QNMs there were twelve possible ordered pairs of modes. Figure El shows a as a function 
of p for four different signal-to- noise ratios. Each a(p) curve is constructed from 10 4 
simulated observations with that amplitude-squared signal-to-noise in each mode. 



if 




1 



(19) 
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Figure 6. False alarm probability a as a function of the probability p appearing in 
equation l|15|l . A false alarm is a misidentification of a QNM pair as arising from 
something other than a general relativistic black hole. 



4-3. False dismissal probability calculation 

Complementary to a, the probability that we incorrectly decide we have observed QNMs 
from something other than a black hole, is the probability that we falsely conclude we 
have observed QNMs from a black hole. This probability is referred to as the false 
dismissal probability and commonly denoted (3. 

The false dismissal probability depends on the detailed character of the source, 
which is not a black hole. Strong gravitational wave sources are compact, with radius R 
on order their mass GM/c 2 and oscillations periods of order GM/c. At the frequencies 
where LISA will have its greatest sensitivity — 10~ 2 -10~ 4 Hz, corresponding to masses 
of order 10 6 -10 8 M Q — we know of no compact sources that are not black holes. For 
the purpose of illustration and to give a sense of the ability of the test described here 
to "discover" new physics, we suppose a population of sources whose frequencies and 
damping times share the same relationship as certain neutron star w-modes calculated 
in j2Sj- Referring to table 1, col. 1, lines 3, 5] we consider observations consisting 
of two modes 

Mu x = 0.471 + 0.056i, 

Mu 2 = 0.654 + 0.164*, (20) 
where M is drawn from the distribution given in equation (fTHJ) . In exactly the same way 
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Figure 7. False dismissal probability as a function of false alarm probability /3(a). The 
false dismissal probability depends on the non-black hole QNM spectrum, which we 
have taken to have the same ratio of frequencies and relationship between frequencies 
and damping times as neutron star w-modes. 



that we used simulations in section 14.21 to determine a as a function of p we calculate 
from these simulations (3 as a function of p. Together a(p) and (3(p) determine (3(a), 
which we show in figure [7| A measure of the effectiveness of the test is the degree to 
which the curves for different signal-to-noise fall below the (3 = 1 — a diagonal. (A 
"test" that randomly picked a fraction a of observations as not black holes would have 
(3 = 1 — a. Any "test" that can do better than randomly choosing in this way will have 
a (3(a) curve that falls below this diagonal.) As expected the test also does better with 
stronger signals. Consider a false alarm threshold of 1%. Then for observations with 
S/N equal to 10 we have a better-than-40% chance of distinguishing NBH sources from 
BH sources. This climbs to better-than-90% chance for observations with S/N 100. 



5. Potential for application 

We have shown that, given at least two QNM signals, from the same source and with 
sufficiently large signal-to-noise, we can cleanly distinguish black holes from other 
astrophysical sources. Following |8., and the discussion above we assume that an 
amplitude-squared signal-to-noise ratio of 100 for the weaker mode is sufficiently large. 
In this section we show that LISA detector observations should provide many examples 
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of such observations, permitting the use of this test to uniquely identify the existence 
of general relativistic black holes. 

Focus attention on each individual QNM. The signal strength, characterized by the 
signal-to-noise ratio at the detector, depends on the energy radiated in the mode, the 
radiation pattern associated with the mode, and the relative orientation of the detector 
and the source. We can average over these latter angles to obtain the mean-square 
signal-to-noise associated with the nim mode as a function of the mode energy j^Hl eq. 
2.30]: 

2 2(1 + z) 2 [°° 1 dE 

where z and D(z) are, respectively, the redshift and the luminosity distance to the 
source. 

For the QNM, we take 

hnim = Antra exp ( - ^"^M s i n {2l\f n n m t) (22) 



Qnim J 
Where Q n lm = ^fnimTnlm- 

Note that Q n e m , which is an observable property of a QNM, is independent of source 
redshift, while the observed f n e m and r n £ m depend on redshift. 

The ringdown energy spectrum of the nlm mode is taken from eq (3.18) of [3*U] : 

dE e £n£m Qnlm -M f 



df e F nlm (AQ 2 n£m + 1) 7T 3 r 2 L [(/ - fnim) 2 + (2ttt 



,-212 



+ 



[(/ + /n. m ) 2 + (27rr)-2p_ (23) 

where the mode amplitude A nim has been replaced with the fraction e n £ m of the mass 
radiated in that mode, defined by 

Using this spectrum in the formula above, and assuming S n (f) is constant over the signal 
band we can integrate over frequencies and invert the result to obtain an approximate 
range over which we can observe a mode n£m with signal-to- noise > p 2 £ m : 

T)(A 2 = 8 Qnim (! + z) 3 M 3 e nlm G 3 
U 57^4^ + 1 F 2 ntm S nP l £m cT [b) 

Given a threshold p 2 nim , black holes radiating a fraction e n £ m of their rest energy in 
mode n£m are observable with a redshift z satisfying eq (|23j) . To use this relation, we 
must make some reasonable assumptions about the observed modes: 

Numerical simulations (eg., |41j ) suggest that energy emitted in QNMs during 
ringdown can be as high as 3% of the rest-mass energy of the hole. For equal-mass 
black hole mergers, the simulations suggest that the i = 2 modes will be by far the 
strongest, with total emitted energies greater than I = 4 modes by as much as three 
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orders of magnitude (see 40j). Here we assume that the weaker mode of a QNM pair 
carries away a fraction 10~ 5 of the black hole mass. 

For LISA, the noise power spectral density is expected to be least in the 10~ 3 < 
/ < Hz band, where it is estimated to be 5 x 10~ 45 Hz -1 . We take a frequency 
/ = 10~ 2 Hz here, with this minimum noise power spectral density as S n . 

The relationship between luminosity distance and redshift we take to be given by eqs 
(23,25) of [37]. This depends on cosmological parameters H (the Hubble parameter), 
Qm, &a and Qk', we use values from a recent review [HE]- 

Finally we can estimate the redshift within which we can expect that supermassive 
black hole mergers will be visible. In addition to the above assumptions, we require a 
signal-to-noise in the weaker mode of at least p 2 = 100. Then, assuming both modes 
are visible in the LISA band, LISA will observe the "average" merger within a redshift 
of ~ 475 (for extremal-spin Kerr) or ~ 320 (for zero-spin Kerr), that is, the entire 
observable universe. 

The rate of merger of black holes of the appropriate mass depends on redshift, both 
because of evolution, and because the observed frequency scales as [(1 + z)M]~ l . The 

estimates of the event rate out to, say, z — 1, then range from 0.1 100/yr (see 

|42j . |4*3~] for recent reviews). 

We conclude that, if any mergers are observable, there is good reason to believe 
they will result in strong QNM excitation and several modes will be observable and 
separately identifiable. 

6. Conclusion 

We have described a qualitatively new test of the existence of general relativistic black 
holes, based on the gravitational radiation they emit when they are formed or when they 
are impulsively excited, e.g., through a merger event. Radiation from an impulsively 
excited black hole, such as might arise in the course of a non-spherical black hole 
formation event or the coalescence of a black hole with another black hole or compact 
object, has a component that consists of a sum of damped sinusoids. This signature 
is, characteristic of the radiation from any impulsively excited, damped source. For 
any given mode, the scale of the frequency and damping time measures the black hole 
mass and angular momentum. Similarly, the relationship of the different modes to each 
other — i.e., the spectrum — is unique to black holes. We have described here how 
this relationship can be used to test the proposition that observed gravitational waves, 
characteristic of an impulsively excited, damped source, in fact originate from a general 
relativistic black hole. Such a test can be characterized in at least two different ways: 
as a definitive "proof" that a black hole has been observed, or as a test of the so-called 
"no-hair" theorem of general relativity. 

To demonstrate the effectiveness of this test we have evaluated numerically the 
probability that the test will mistakenly fail to identify an actual black hole. By 
introducing a hypothetical gravitational wave source whose characteristic frequencies 
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and damping times are similar to those of neutron star w-mode [221 we have also 
evaluated numerically the probability that the test will incorrectly identify w-mode 
oscillations of a neutron star, or any object whose spectrum is similar to a black hole. 
Together these results demonstrate that for sources with the signal-to-noise expected 
of, for example, massive black hole coalescences detected by LISA, the test proposed 
here can cleanly discriminate black hole sources. Finally, we have shown that LISA can 
be expected to observe signals of this kind and strength throughout the universe with 
a rate that may be as great as 100/y. 

This method can be used to measure mass and angular momentum of a black hole. 
Using gravitational waves to measure mass and angular momentum is an idea that has 
been around for some time [JJ |S] . In these previous works it was assumed that the mode 
observed was of a known order (e.g., the mode with the longest damping time, or the 
lowest order, etc.). With the observation of two or more modes the requirement that 
a single mass and angular momentum explain the complete set likely permits the mass 
and angular momentum to be determined uniquely. 

The field of gravitational-wave detection is new. The current generation of ground- 
and space-based gravitational-wave detectors is opening a new frontier of physics: 
gravitational-wave phenomenology, or the use of gravitational wave observations to learn 
about the physics of gravitational-wave sources and gravity itself. We are only just 
beginning learn how to exploit the opportunities it is creating for us. As gravitational- 
wave observations mature, we can expect more and greater recognition of their utility as 
probes of the character of relativistic gravity. The opening of this new frontier promises 
to be an exciting and revealing one for the physics of gravity. 
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